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Abstract 

In  this  paper  we  consider  B-wavelets  of  order  2,  i.e.  piecewise  linear  spline  prewavelets  of 
smallest  support,  over  nonuniform  knot  sequences.  We  discuss  an  example  showing  that 
for  1 < p < oo,  there  is  no  absolute  L;, -stability  for  these  B-wavelets.  This  means  that 
regardless  what  specific  scaling  of  the  B-wavelets  is  chosen,  the  corresponding  stability 
constants  cannot  be  made  independent  of  the  knot  sequences  involved. 


1 Introduction 

Polynomial  splines  are  fundamental  tools  in  numerous  branches  of  applied  mathematics, 
and  for  spline  spaces  defined  over  a given  knot  sequence,  the  basis  of  choice  is  provided  by 
B-splines,  which  possess  a lot  of  attractive  properties  for  numerical  computations.  One 
of  these  important  properties  of  B-splines  is  their  absolute  stability.  Given  a B-spline 
basis  {Bi}ieI  of  polynomial  order  d over  a valid  knot  sequence  t,  a classical  result  by  de 
Boor  [1]  states  that  properly  normalized  B-splines  are  stable  in  the  sense  that  for  each 
set  {bi}if:T  of  real  coefficients  it  holds  that 

OMIp< 

Here  ||  • ||p  denotes  the  standard  integral  and  discrete  p-norms  for  1 < p < oo,  respectively, 
and  the  normalizing  factor  Si  for  each  B-spline  is  the  length  of  its  support  divided  by 
the  order  d.  The  important  point  is  that  the  positive  constant  C<j  is  dependent  on  the 
order  d alone,  and  not  in  any  way  on  the  underlying  knot  sequence  t. 

Since  nested  knot  sequences  give  rise  to  nested  spline  spaces,  spline  functions  have  also 
become  a focus  of  attention  within  the  theory  of  wavelets  and  multiresolution  analysis, 
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starting  with  cardinal  spline  wavelets  on  infinite  equally  spaced  and  uniformly  refined 
knot  sequences,  for  which  Fourier  transform  techniques  are  available,  see  [3]  and  the 
references  therein. 

The  study  of  spline  wavelets  on  bounded  intervals,  for  arbitrary  knot  sequences  and 
nonuniform  refinement  began  with  the  papers  [4],  [5]  and  [2],  respectively.  The  construc- 
tion of  so-called  minimally  supported  B- wavelets  for  a given  spline  order  d and  two  nested 
knot  sequences  to  provide  a basis  of  the  relative  orthogonal  complement  (wavelet)  space 
is  described  in  detail  in  [6].  This  means  that  given  the  coarse  and  fine  knot  sequence, 
there  exist  explicit  algorithms  to  determine  the  supports  of  the  B-wavelet  functions,  the 
so-called  minimal  intervals,  and  also  to  compute  the  corresponding  wavelet  functions, 
though  only  up  to  a normalization  constant. 

One  open  problem,  however,  is  how  to  fix  the  normalization  factor  for  each  B-wavelet 
function  to  achieve  best  possible  stability  for  the  whole  B-wavelet  basis.  We  provide  an 
example  for  the  case  of  piecewise  linear  wavelets,  i.e.  polynomial  order  2,  that  shows 
that  for  1 < p < oo  there  is  no  absolute  stability  of  B-wavelets,  meaning  that  there  is  no 
choice  of  normalization  that  provides  absolute  stability  constants  which  are  completely 
independent  of  the  underlying  knot  sequences.  ip-stability  estimates  involving  a quantity 
dependent  on  the  knot  sequences  for  1 < p < oo  and  showing  absolute  stability  for  p = 1 
are  given  in  [7]. 

2 Piecewise  linear  B-wavelets 

The  theory  of  B-wavelets  [6]  covers  general  cases  of  knot  refinement,  such  as  situations 
where  several  or  no  knots  at  all  are  inserted  into  an  old  knot  interval,  or  where  the 
multiplicity  of  an  existing  knot  is  increased.  For  our  purposes,  however,  it  is  sufficient  to 
consider  what  one  might  call  the  standard  setting,  where  all  knots  are  simple  except  at 
the  interval  endpoints,  which  we  can  count  as  double  knots,  and  where  exactly  one  new 
knot  is  inserted  strictly  between  two  old  ones. 

Our  notations  are  as  follows  for  the  closed  interval  [0, 1] . We  have  a coarse  knot 
sequence  with  n—  1 interior  knots,  namely 

r : 0 = To  < Ti  < • • • < r„  = 1. 

Strictly  between  each  pair  of  coarse  knots  t<_  i and  Ti  we  insert  a new  knot  Si  at  an 
arbitrary  location,  i.e. 

n- 1 < Si  <Ti  for  each  i = 1, . . . , n. 

Thus  we  have  a sequence  s of  new  knots 

s : 0 < Si  < • • • < s„  < 1. 

The  fine  knot  sequence  t = r U s,  when  ordered  appropriately,  is  given  as 

t : 0 = to  < tl  < • • • < t2n  — 1) 

where  the  even  numbered  knots  in  t correspond  to  old  knots  in  r,  while  the  odd  numbered 
knots  represent  the  newly  inserted  knots  from  s.  To  account  for  the  boundary,  we  treat 
the  interval  endpoints  as  double  knots  by  setting  r_i  = t_i  = 0 and  rn+\  = t2n+:  = 1- 
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For  our  investigations  it  is  necessary  to  introduce  also  some  notation  related  to  the 
knot  spacings.  We  set 

di  — ti+ 1 - ti  for  t = 0, . , 2n  — 1,  and  <5;  = U+\  - f;_i  for  i = 0, . . , , 2 n, 


which  means  So  = do  = h - to  and  62 n = d2n-i  = t2n  — hn-i  at  the  boundary.  Thus  St  is 
the  distance  between  two  consecutive  old  knots  if  i is  odd,  and  between  two  consecutive 
new  knots  if  i is  even  (and  not  at  the  boundary). 

We  also  introduce  the  index  sets 

0 - {1, 3, . . . , 2n  - 1}  and  fi0  = {3, 5, . . . , 2n  - 3}  . 

The  piecewise  linear  functions  on  the  knot  sequences  r Ct  form  nested  linear  spaces 
Vo  C V\  of  dimensions  n+ 1 and  2n+l,  respectively.  The  corresponding  piecewise  linear 
B-splines  forming  a basis  of  these  spaces  are  simple  hat  functions.  We  denote  them  as  <pj 
and  7 i for  r and  t,  respectively,  where  with  the  necessary  adjustments  at  the  endpoints, 


(x  ~ Tj- 1)  /S2j- 1 if  X £ [ Tj-iTj ] 


<Pj  (x)  = < 

, (u+i- 

■ £)  /f>2j+l 

if  X e [Tj,Tj+ 1]  for  j = 0, . . . , n, 

(2.1) 

1 

L 

0 

otherwise 

( 

-1)  /di- 1 

if  x £ [U-iti] 

7 i (x)  = < 

j (ti+ 1 

-x)  /di 

if  x £ [tj.tj+ij  for  i = 0, ... , 2 n. 

(2.2) 

■ '' 

l 

0 

otherwise 

Using  for  any  two  functions  /,  g € Vj  the  standard  inner  product 

(f,9)=[  f(t)g(t)dt: 

Jo 

we  can  write 


Vi  = Vo  © W, 

where  W is  the  relative  orthogonal  complement  of  Vo  in  Vi , and  0 denotes  orthogonal 
summation.  The  dimension  of  W is  n,  so  that  there  is  a basis  function  ipk  for  every  index 
k £ O,  or  in  other  words  for  each  newly  inserted  knot  s*. 

Nonzero  functions  ipk  € W with  minimal  support  are  called  B-wavelets.  The  general 
theory  for  B-wavelets  developed  in  [6]  establishes  in  this  special  case  that  there  are  n 
different  piecewise  linear  B-wavelets  which  form  a basis  of  the  wavelet  space  W . Each 
such  B-wavelet  is  uniquely  determined  up  to  a constant  multiple.  There  are  two  boundary 
B-wavelets  ipi  and  fan- 1 and  n — 2 interior  B-wavelets  ipk  for  k € fio>  which  we  will 
consider  first.  Each  interior  B-wavelet  has  support  [^-3,^+3],  so  that 

k+ 2 

fa  (z)  = dtli  (*)  for  x € [0,  1] 

i—k—2 

with  the  coefficients  determined  by  ipk  € W,  or  in  other  words 

{ipk,<pj)  = 0 for  j = 0, 
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For  the  boundary  wavelets  ipi  and  ip2n-i  we  have  to  make  some  minor  modifications. 
Their  supports  are  [tot^]  and  [t2n-4t2n],  respectively,  so  that 

3 2 n 

ipi  (x)  = ^ qjji  ( x ) and  if)2n-i  ( x ) = q?n~lli  {x)  for  x e [0, 1] . 

i— 0 i=2n— 3 

In  the  paper  [7]  the  values  of  all  B- wavelet  coefficients  qf  are  given  explicitly  in  terms 
of  the  knot  locations  for  the  standard  setting  described  here.  In  the  same  paper  estimates 
for  the  coefficients  are  used  to  derive  ip-stability  estimates  for  these  B- wavelets. 

3 Stability  of  B-wavelets 

Our  aim  in  this  paper  is  to  establish 

Theorem  3.1  Given  the  B-wavelet  basis  {tpk}ken  > then  for  1 < p < 00,  there  are  no 
sets  of  weights  ak,p>  k € Cl,  such  that 

Ki\\c\\p<  £c*afclPife  <^2||c||„  " (3.1) 

keo.  p 

holds  for  any  wavelet  coefficients  (c*, C3, . . . , C2n-i)  and  with  absolute  constants  K\  > 0 
and  K2  > 0,  which  are  completely  independent  from  the  choice  of  knot  sequences  r and  s. 

Due  to  the  finite  dimension  of  W,  it  is  clear  that  stability  constants  Kt  and  K2  exist, 
as  any  two  norms  on  W are  equivalent.  The  pertinent  question  is  how  the  weights  could 
be  chosen  to  achieve  that  the  constants  are  actually  independent  of  the  dimension,  the  p- 
norm  and,  if  possible,  the  choice  of  new  knots  s.  We  will  prove  the  assertion  by  assuming 
that  the  estimate  (3.1)  holds  with  constants  independent  of  the  knot  sequences.  Then 
the  following  special  case  serves  as  a counterexample  to  this  assertion. 

The  old  knot  sequence  r consists  of  the  equally  spaced  points: 

r0  = 0,Ti  = 1/3, r2  = 2/3,  t3  = 1. 

We  want  to  investigate  what  happens  if  two  newly  inserted  points  are  positioned  ever 
more  closely,  so  we  introduce  the  new  knots  as 

si  = 1/3  — e,S2  = 1/3  + 77,  S3  = 5/6,  for  0 < e,  q < 1/3, 

in  order  to  find  out  what  happens  if  both  e — > 0+  and  77  — > 0+. 

Thus  the  fine  knot  sequence  t is 

t0  = 0 , ti  = 1/3  -e,t2  = 1/3,  h = 1/3  + 77,  tA  = 2/3,  t5  = 5/6,  f6  = 1. 

The  fine  interval  lengths  are 

do  — 1/3  - s,  d\  = e,  d2  = 77,  ds  = 1/3  - 77,  d4  = 1/6,  do  = 1/6, 

l 

di  = S3  — S5  = 1/3,  and  <f0  = 1/3  - s,S2  = £ + 77, 64  = 1/2  - 77,  <f6  = 1/6. 


while 
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In  this  setting  any  wavelet 


t= 0 


must  be  orthogonal  to  the  coarse  hat  functions  <p0,  ■ ■ ■ , ip3.  This  actually  means  that  the 
column  vector  q of  coefficients  qt  must  satisfy  the  matrix  equation 

Aq  - 0,  (3.2) 

where  the  entries  of  A are  the  inner  products  of  the  coarse  and  fine  hat  functions,  i.e. 

Oj,*  = (¥b>7i),  for  j = 0, . . . , 3,  i = 0, . . . , 6. 

Direct  computations  using  (2.1)  and  (2.2)  yield  as  the  only  nonzero  entries 


0-0,0 

0-1,0 

0-1,2 


loll  11  1 o 

"2£  _6£+9’  ao’1-6£+18’  a°’2  ~ 2£  ’ 

loll  11 

2£  “3£+l8’  ai-1  = -g£  + 

1 2 1 1 2 1 
"2£  +2£~2V  +2n' 

1 


9’ 


1 


0l-3  = "6^  + i aM  = ^2-^+i8 


1 2 

02,2  = 2^ 


a2'3=67?+18’ 


02A  = -xr?  -T»J  + 


13 


6 ' 72’ 


02,5 

03,4 


02,6  = 


1_ 

12’ 

1 1 

72’  03,5  ~ 12’ 


72’ 


03,6  = 


72 


We  now  investigate  the  B-wavelets  ipi  and  i/;3  in  detail,  corresponding  to  Si  = t\  and 
S2  = t's.  Specializing  the  results  from  [7]  then  yields  all  necessary  B-wavelet  coefficients 
for  this  setting  up  to  a scaling  factor.  Note,  however,  that  it  is  straightforward  to  check 
that  the  corresponding  coefficient  vectors  satisfy  the  matrix  equation  (3.2). 

The  coefficients  of  the  boundary  wavelet  ipi  are 

3 

1 — 3e’ 

6_  9eri 

£ + T)  + 6 £Tj  ’ 

1 + 3r? 
e + r]  + 6er]  ’ 

9r f 

e + T]  + 6er 7’ 


q\ 

qI 
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while  the  ones  for  the  interior  B-wavelet  ip 3 are 


4 = 


4 


9e2 


q3  — 3 + 


4a  = 


£ + T]  + 6 £7?  ’ 

l + 3s 

£ + r)  + 6 er)  ’ 

3(e  + V) 


+ 


9(1  ~ 277) 
2(£  + ??  + 6£?})  ' 2(5-12)7)’ 
9 


«4 


5-1277’ 

3 


2 (5  - 1277) ' 

We  first  provide  estimates  for  the  p-norms  of  these  B-wavelets. 
Proposition  3.2  For  small  enough  e and  77,  it  holds  for  1 < p < 00  that 

16 


IIV’i 


> 

'p  - 45 


(ff<-^- 


16/1\ 


i/p 


llfcllp  > (*  + »?) 

Proof:  For  all  0 < e,  77  < 1/3  we  find  that 


l/p-i 


|o5l  > (»+7r1,<M/,!‘,'^,^) 
- 1^+^' 


and,  similarly, 


92 1 > g (e  + v) 


-1 


Note  that  instead  of  8/9  we  may  write  1 — <7  for  any  a > 0 if  e and  77  are  small  enough 
or  even  1 if  e = 77. 

In  the  process  e,  77  — *•  0+  all  other  coefficients  q}  and  qf  have  finite  limits.  This  means 
that  for  small  enough  e + 77 


8 


max \q}\  = \q\  | > - (e  + 77)  , 


IIV>3 


max  | = |<?f|  > - (e  + 77)  1 . 


The  absolute  stability  of  piecewise  linear  B-splines  (1.1)  yields  with  C2  > 5/2  (see  [1]) 
and  S2  = £ + 77 


llVhllp  = 


i= 0 


i/p 


5 \ 2 


(?0^0/P>  ?Ml1/P>  92^2/P5  <l35l/P) 
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Analogously  we  get 


ll«p=  X>b.  >?(i)  kll^/p a 3§(l)  (E' 

i—  1 n 


to  complete  the  proof. 


Proof  of  Theorem  3.1:  Let  us  now  assume  that  with  some  scaling  factor  B-wa,velets 
are  absolutely  stable  in  p-norm  for  1 < p < 00,  i.e.  there  exist  weights  a*-,p  so  that 
the  inequalities  (3.1)  hold  with  constants  independent  of  the  specific  choice  of  knot 
sequences.  Choosing  in  the  current  setting  all  coefficients  equal  to  zero  except  for  c\  = 1, 
the  stability  inequality  (3.1)  yields 

IK, Philip  < k2 

or  in  other  words,  using  Proposition  3.2 

|QiJ£ra;sl2,/,ft(£+’,)1"/'  (3-3) 

and  by  a similar  argument 

M£m;£^2'/rK2{c+’,)'~'/r'  (34) 

On  the  other  hand,  the  stability  estimate  (3.1)  yields  for  arbitrary  C\  and  c 3,  while 
setting  all  other  c*,  to  zero,  that 

\\ciai,pijh  +c30i3ipV’3||p  > Ki  (|ci|p  + MP)1/P  > /Limax(|d| , |c3|). 

Let  us  choose  specifically 

Ci  = a3iP  and  c3  = —a iiP, 

which  results  in 

|ai,pa3,p|  HV'i  - 1P3 ||p  > ATimax  (|qi>p|  , |a3iP|) , 
leading  with  (3.3)  and  (3.4)  to 


(£  + r?)1-1/p||V’i-V'3||p>  Q)1 


On  the  other  hand  we  derive  from  the  absolute  stability  of  linear  B-splines  (1.1)  that 
= 1 1 <?o To  + (9}  - 9?)  7i  + (92  - 9f)  72  + (93  - 93)  73  - 9474  - qha\\P 

<||  {iX'r,  (d  - «?)  «!/p,  (4  - 4)  4/p.  (95  - « 4/p.  9K/p.9$/p)  ||p 


An  Lp-stability  example  for  piecewise  linear  B- Wavelets 
< 6max(|gi|<So/p,  |«i  — 9i| 

All  the  terms 

kol  <*d/p,  \q\  - «?|  <5i1/p5  ks  - qI\ 4/p>  k!|  <5l/p>  loti  4/p 

are  in  fact  bounded  from  above  for  e + rj  — * 0+ , so  that  the  expression 

(e  + f?)1_1/P|9oK/p 

and  the  other  such  terms  tend  to  zero. 

Since  \q\  — q% I = 3 |e  — rj\  / (e  + rj  + 6 erj)  we  obtain  for  the  only  remaining  term 

(£  ^ - ^1  1 + 6 ev/^+  ,)  - 3 k ~ ' 

which  goes  to  zero  as  well  for  e + r]  — * 0+.  As  a consequence 

lim  (£  + t?)1-1/p  ||-0i  — V»3||p  = 0, 

£+7}— >0+  y 

which  contradicts  (3.5).  □ 

Remark  3.3  Although  we  have  chosen  an  example  with  one  boundary  and  one  interior 
B-wavelet,  let  us  remark  that  the  lack  of  absolute  stability  is  in  no  way  due  to  a boundary 
effect.  A completely  analogous  reasoning  is  possible  if  one  chooses  knot  sequences  with 
more  interior  knots  and  studies  the  behaviour  for  two  interior  B-wavelets  once  two  new 
knots  coalesce.  Similarly  just  two  boundary  B-wavelets  could  be  used  on  an  even  shorter 
knot  sequence,  where  there  are  no  interior  B-wavelets  at  all. 
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